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Abstract 

Using a large sample of D° — > K~fj, + i/ and D° — > ir~ fi + is decays collected by the 
FOCUS photoproduction experiment at Fermilab, we present new measurements of 
the q 2 dependence for the f+{q 2 ) form factor. These measured f+(q 2 ) form factors 
are fit to common parameterizations such as the pole dominance form and compared 
to recent unquenched Lattice QCD calculations. We find ?7ipoic = 1.93 ± 0.05 ± 
0.03 GeV/c 2 for D° -> K'fi+u and m po i c = 1.91±g;?g ± 0.07 GeV/c 2 for D° -» 

ir-fi+v and f { _ K \o)/f{ K) (0) = -1.7±U ± 0.3. 
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1 Introduction 



In this paper, we provide a new non-parametric measurement of the q 2 evo- 
lution for the f+(q 2 ) form factor describing the pseudoscalar decay D° — > 
K~fi + v. The measurement is presented in a form that is convenient for para- 
metric and non-parametric comparisons to other experiments and theoretical 
predictions. Our q 2 evolution is compared to the lattice gauge calculations in 
[2], and we show that fits to the q 2 evolution agree with traditional parametric 
analyses of the data and results from other experiments. It is important [1] to 
make incisive tests of unquenched Lattice Gauge calculations in semileptonic 
decay processes such as D° — > K~fi + v and D° — > 7T~fi + v in order to ulti- 
mately reduce the substantial systematic errors on the CKM matrix in charm 
and related beauty processes. 

Two form factors describe the matrix element for such decays according to 
Eq. 1 



M = G F V CS f + (q 2 ) (P D + P K ) a + f^q 2 ) (P D - P K ) a utf{l - 7 K (1) 



These lead to a differential width of the form given by Eq. 2 where Pk is 
the kaon momentum in the D° rest frame and all f-{q 2 ) contributions are 
multiplied by the square of the muon mass. 1 
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In Eq. 2, W = (m 2 D + m\ - mj) /(2m D ), F = W - E K +ml/(2m D ) and P K 
, Ek are the momenta and energy of the kaon in the D° rest frame. Assuming 
f-(q 2 )/ f+(q 2 ) is of the order of unity as expected, the corrections due to f-(q 2 ) 
are fairly small and, apart from the low q 2 region, 
dY/dq 2 



G 2 F 



V cq \ 2 P|- |/ + (g 2 )| 2 /(247r 3 ) is an excellent approximation. 



This paper provides new measurements of the f+{q 2 ) form factors for D° — > 
K-fi+v and D° -> t^-^v and of the ratio /_(0)//+(0) for D° -> K~^v 2 
Our emphasis in this paper is on the shape of the q 2 dependence rather than 
on its absolute normalization. As a means of comparing our result to different 
parameterizations commonly used in the literature, we will fit our measure- 
ments of f+(q 2 ) to two different parameterizations in use: the pole form given 



1 This form was obtained using the basic formulae in [3] . 

2 Throughout this paper, we will assume that f-(q 2 )/ f+{q 2 ) is essentially indepen- 

-.2 



dent of q 
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by Eq. 3 and the modified pole form given by Eq. 4. 3 



/ + (o) 



(3) 



1 - q 2 /m 2 



pole 



U(0) 



(4) 



(1 — q 2 /m 2 D *) (1 — a q 2 /m 2 D *) 



Throughout this paper, unless explicitly stated otherwise, the charge conjugate 
is also implied when a decay mode of a specific charge is stated. 



2 Experimental and analysis details 

The data for this paper were collected in the photoproduction experiment 
FOCUS during the Fermilab 1996-1997 fixed-target run. In FOCUS, a forward 
multi-particle spectrometer is used to measure the interactions of high energy 
photons on a segmented BeO target. The FOCUS detector is a large aperture, 
fixed-target spectrometer with excellent vertexing and particle identification. 
Most of the FOCUS experiment and analysis techniques have been described 
previously [4,5,6]. In this section we describe the cuts used both in the non- 
parametric analysis described in Section 3 as well as the parametric analysis 
described in Section 4. 

The non-parametric part of this analysis is based on a sample of ~ 13, 000 
decays of the form D* + — > D°ir + , where D° — > K~fi + v. To isolate the D° — > 
K~ n + v topology, we required that the muon, and kaon tracks appeared in a 
secondary vertex with a confidence level exceeding 1%. In order to suppress 
backgrounds from higher multiplicity charm decays, we isolated the K~fi + 
vertex from other tracks (not including tracks from the primary vertex) by 
requiring that the maximum confidence level for another track to form a vertex 
with the candidate be less than 1%. The D* decay pion was required to lie in 
the primary vertex. 

The muon track, when extrapolated to the shielded muon arrays, was required 
to match muon hits with a confidence level exceeding 1% and all other tracks 
were required to have confidence level less than 1%. The muon candidate was 
allowed to have at most one missing hit in the 6 planes comprising our inner 
muon system and a momentum exceeding 10 GeV/c. In order to suppress 

3 In this form, the parameter a gives the deviation of f+(q 2 ) from spectroscopic 



pole dominance where m£)* = mp* = 2.112 GeV/c 2 for D° — > K [i + v and mp* = 
m D+ » = 2.010 GeV/c 2 for D° -» n'fi + u 



4 



muons from pions and kaons decaying within our apparatus, we required that 
each muon candidate had a confidence level exceeding 1% to the hypothesis 
that it had a consistent trajectory through our two analysis magnets. The 
kaon was required to have a Cerenkov light pattern more consistent with that 
of a kaon than that of a pion by 1 unit of log likelihood [6]. Non-charm and 
random combinatorial backgrounds were reduced by requiring a detachment 
between the vertex containing the D° —>■ K~ /i + z/ and the primary production 
vertex of at least 5 standard deviations. 

Possible background from D° — > K~tt + , where a pion is misidentified as a 
muon, was reduced by requiring the reconstructed K~{i + mass to be less than 
1.812 GeV/c 2 . Finally we put a cut on the confidence level (CL closure ) that the 
event was consistent with the hypothesis D° — > K~fi + u that will be described 
below. 




0.14 0.15 0.16 0.17 0.18 0.19 
M(K~/j, + v tt + )-U(K~/j, + v) (GeV/c 2 ) 



Fig. 1. The D* — D mass difference distribution for events satisfying our signal 
selection cuts. The points with error bars are for wrong-sign subtracted data. The 
histogram shows the mass difference distribution obtained in data for the wrong sign 
events. The solid curve shows the wrong-sign subtracted 5m distribution obtained in 
Monte Carlo; while the dashed histogram shows the 5m distribution of wrong-sign 
Monte Carlo events. There is an excess of 12,840 opposite charge combinations over 
same charge combinations where the D* — D mass difference was less than 0.160 
GeV/c 2 . 

The 5m = m(K~[i + v 7r + ) — m(K~ fi + v) distribution for our tagged D° — > 
K~fi + u candidates is shown in Figure 1. Several of the distributions shown 
in Figure 1 are wrong-sign subtracted meaning that combinations where the 
D* decay pion have the same charge as the kaon are subtracted from those 
where the decay pion and kaon have the opposite charge. Figure 1 was cre- 
ated using our standard [7] line-of-flight neutrino closure technique. Briefly, 
the standard neutrino closure method assumes the reconstructed D momen- 
tum vector points along the displacement between the secondary and primary 
vertex. This leaves a two-fold ambiguity on the neutrino momentum. For Fig- 
ure 1, we use the neutrino momentum that resulted in the lower 5m. Figure 1 
illustrates the necessity of making a wrong-sign subtraction, since the wrong- 
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sign fraction in the data is much larger than predicted by our charm Monte 
Carlo, indicating the presence of a non-negligible non-charm background in 
the data. We believe that this non-charm background is right-sign, wrong-sign 
symmetric. 

We next describe the cuts used for the two-dimensional fit analysis described 
in Section 4. One of the principal motivations for this analysis, was to compare 
the decay widths for D° — > it~ n + v and D° — > K~{i + v and extract the ratio 
/+^(0)//+^(0). As such, the analysis cuts used for the two-dimensional anal- 
ysis are somewhat different than those previously described for the deconvo- 
lution analysis in order to reduce systematic uncertainty on the ratio. Several 
additional cuts on the muon candidate were applied to remove contamination 
from electrons. The kaon, in D° — > K~fi + u, was required to have a Cerenkov 
light pattern more consistent with that for a kaon than with that for a pion 
by 3 units of log likelihood. The pion track, in D° — > n~fi, + is, was required to 
have a Cerenkov light pattern more consistent with that of a pion than that 
of a kaon by 3 units of log likelihood. For both the pion from D° — > 7r~/i + z/ 
decay as well as the D* + decay pion, we further required that no other hypoth- 
esis was favored over the pion hypothesis by more than 6 units of likelihood. 
The pion in D° — > n~ fi + u was required to have a momentum greater than 14 
GeV/c , and the D* decay pion was required to have a momentum greater 
than 2.5 GeV/c. A D* — D mass difference cut of 5m < 0.154 GeV/c 2 was 
applied. Finally the hadron-muon mass was required to exceed 1 GeV/c 2 . 

To improve the q 2 resolution for both the parametric and non-parametric 
analyses, we developed an alternative neutrino closure technique that we will 
call D* cone closure. We require that the K~ fi + v reconstructs to the mass of 
a D° and the K~fi + u n + reconstructs to the mass of D* + . When viewed in the 
K~fi + rest frame, these constraints place the neutrino momentum vector on a 
cone about the D* decay pion where both the neutrino energy and cone half- 
angle are determined from the mass constraints and the well measured K~ , 
and 7r + momentum vectors. We then sample all azimuths for the neutrino 
in this cone, reconstruct the lab frame D° momentum vector, and choose the 
azimuth where the D° is most consistent with pointing to the primary vertex 
based on minimizing a x 2 variable. In order to further reduce backgrounds, 
we required CL closurc > 1% where CL closure is a confidence level based this 
minimal \ 2 ■ 

Averaged over all detected D° — > K~fi + v events, the Monte Carlo predicted a 
rather non-Gaussian q 2 resolution with an r.m.s. width of 0.22 GeV 2 /c 2 using 
the D* cone closure technique. The correlation between the generated and 
reconstructed q 2 is illustrated in Figure 2. 

It was important to test the fidelity of the simulation with respect to the 
reproducibility of the q 2 resolution. To do this, we studied tagged, fully- 
reconstructed D° — > K~7t + 7t + 7t + decays from D* + — > D°tt + where, as a test, 
one of the D° decay pions was reconstructed using our neutrino cone closure 
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Fig. 2. A scatter plot of q 2 reconstructed by the cone closure technique and the 
true, generated q 2 in Monte Carlo events. We show the line where q 2 generated 
equals q 2 reconstructed. 

technique. We then reconstructed the q 2 using the neutrino closure and com- 
pared it to a precisely reconstructed q 2 obtained from the magnetically recon- 
structed "neutrino" pion. The difference between these two q 2 values provided 
a resolution distribution obtained from data that could then be compared to 
the same resolution distribution obtained using tagged D° — > fT~7r + 7r + 7r + in 
our Monte Carlo. The Monte Carlo resolution distribution was a good match 
to the observed resolution distribution. 



3 Non-parametric analysis of /+ \q 2 ) 

We begin with a description of the method used to correct for the effects of 
acceptance and q 2 resolution. We will call this our deconvolution technique. 
The goal of the deconvolution is to produce a set of f 2 (q 2 ) that represents 
measured f'^ < \q 2 ) values - each averaged over the (narrow) width of the q 2 
reporting bins. Under the assumption that f-(q 2 )/ f+(q 2 ) is on the order of 
unity, Eq. 2 implies that the number of events expected in a given q 2 bin is 
proportional to \f+(q 2 )\ 2 - Our Monte Carlo is used to determine the fraction 
of events reconstructed in a given q 2 bin that were generated in another q 2 
bin. This information, along with the f'^ < \q 2 ) distribution used in the original 
generation 4 of the D° — > K~fi + u Monte Carlo sample, was combined to form 
a square matrix that linearly relates a vector of the predicted number of events 
reconstructed in each q 2 bin to a column vector of assumed f1(q 2 ) values. 




(N n /fl(q 2 ) N 12 /f 2 + (q 2 ) ...\ 
N^/PM) N 22 /p + {q 2 ) ... P + {q 2 ) 





4 The sample was generated assuming f-{q 2 )/ f+(q 2 ) = —0.7 
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In Eq. 5, rii are the number of events that are observed in the ith q 2 bin, rj is 
the ratio of the observed to the number of generated events, JV^ is the number 
of Monte Carlo events that were generated in q 2 bin j that reconstruct in q 2 
bin i, f 2 (q 2 ) are the input f 2 (q 2 ) used in the Monte Carlo generation, and 
f+(qf) are the true form factors that describe the data. 

With reference to Eq. 5, the "deconvolved" f+{qf) is then given by the inverse 
of the square matrix times the first column vector that consists of the observed 
number of events reconstructed in each q 2 bin. We will call the inverse of 
this matrix the deconvolution matrix with components D^. In this notation: 
PM) = £7 Aj n r 

We performed the sum f 2 (q 2 ) = J2j D^rij by using a separate, weighted 
histogram for each f+(q 2 ). Each / 2 (g 2 ) is a sum of weights over all events 
where the event weight Dij where % is q 2 reporting bin and j is the reconstructed 
q 2 bin for that event. We perform a wrong-sign subtraction by multiplying the 
deconvolution weight by +1 if the kaon of the event had the opposite 
sign of the D* decay pion and -1 otherwise. Our charm background correction 
is based on a Monte Carlo, which incorporates all known charm decays and 
charm decay mechanisms. The charm background was normalized to the same 
number of D* + — > D°tt + , D° — > K~tt + events observed in the data. We 
subtract known backgrounds by deducting the deconvolution weights for the 
background events predicted by our charm Monte Carlo. Figure 3 compares 
the f+(q 2 ) values obtained with and without the background subtraction. 5 
Figure 3 shows that most of the charm background is expected in the high 
q 2 region, and once the background is subtracted, the data is an excellent fit 
to the pole form. In the range 5m < 0.16 GeV/c 2 , the expected, wrong-sign 
subtracted background yield from our Monte Carlo was found to be 12.6 % of 
the total number of events in this 5m range when using our baseline cuts. 

The deconvolution was obtained by summing the weights of all events with 
5m < 0.16 GeV/c 2 . A ten bin deconvolution matrix was used, with the over- 
flow bin dropped. Given that our bin width, 0.18 GeV 2 /c 2 , is comparable to 
our r.m.s. q 2 resolution, adjacent f+{q 2 ) values have a strong, negative corre- 
lation (typically - 65%) and the error bars are thus significantly inflated over 
naive counting statistics errors. Table 1 summarizes of our non-parametric 
f^\q 2 ) measurements for D° — > K~fi + u along with the correlation matrix. 
The f^iq 2 ) values are normalized such that /+(0) = I. 



5 The covariance between two q 2 reporting bins is given by the sum of the product 
of event weights for the two reporting bins. We take the square root of the f+(qf) 
returned by the fit and make the appropriate adjustment to the variances obtained 
from the diagonal elements of the covariance matrix. 
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0.5 
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0.5 1 1.5 2 

q 2 (GeV 2 /c 2 ) 

Fig. 3. The deconvolved f^\q 2 ) for D° — > K~fi + v events using nine, 0.18 GeV 2 /c 2 
bins. The triangular points are before the known charm backgrounds were sub- 
tracted. The square points are after subtraction for known charm backgrounds. The 
line shows the pole form with m po i c = 1.91 GeV/c 2 - the value obtained from a fit 
to the displayed /^(q 2 ) points. After background subtraction, the confidence level 
of the fit to the pole form is 87%. A fit to the modified pole form produced an a 
parameter of 0.32 with a confidence level of 82%. 



Table 1 

Measurements of f^\q 2 ) for D° — » K^fi^u along with the 9x9 matrix of relative 
correlation coefficients. 
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4 Parameterized f+(q 2 ) forms for D° — > K and D° — > 7r /x + z/ 

In this section we present values of the m po i c and a parameters for the pole 
form (Eq. 3) and modified pole form (Eq. 4) as well as the ratio (0) / (0) . 
For the case of D° — > K~fi + u we have done this directly from a x 2 fit of 
the non-parametric f^\q 2 ) values illustrated in Figure 3 as well as from a 
two-dimensional binned likelihood fit to the q 2 and cos 6^ scatterplot where 
0£ is the angle between the v and the D direction in the \xv rest frame. Be- 
cause of the much larger background contamination in the Cabibbo suppressed 
D° — > 7T~/i + z/, only the two-dimensional binned likelihood fit was employed to 
extract f^\q 2 ) parameters for this mode. We begin with a discussion of the 
results from the two-dimensional fit. 

Information on parameterized f^\q 2 ) and /i^(0)//+^(0) is obtained by us- 
ing a weighting technique that is similar to that described in [7]. We use a 
binned version of the fitting technique developed by the E691 collaboration [8] 
for fitting decay intensities where the kinematic variables that rely on recon- 
structed neutrino kinematics are poorly measured. The observed number of 
events in each g 2 -cos 9g bin is compared to a prediction based on signal inten- 
sity as well as background contributions. 

The signal component is constructed from a weighted Monte Carlo. The signal 
Monte Carlo was initially generated using nominal values for /_(0)// + (0) and 
mpole- Both the generated as well as the reconstructed kinematic variables 
were stored for each event. The signal prediction for a given fit iteration is 
then computed by weighting each event within a given reconstructed kinematic 
bin by the intensity evaluated using the generated kinematic variables for the 
current set of fit parameters divided by the generated intensity. 

A variety of possible backgrounds were included for our two processes. These 
included general charm background based on our charm Monte Carlo as well 
as specific backgrounds that create peaks in the 5m distribution. For the 
case of D° — > n~fi, + is, the specific backgrounds included K~fi + v, K~7r fi + v, 
K 7r~ /i + z/ and p^fiu ; while for D° — > K~ /x + z/ this included K~ "ttVz/. In all 
cases, the shape of the backgrounds were determined from our Monte Carlo 
that incorporated known decay intensities. The branching ratio of each specific 
background relative to the two signal processes were allowed to float, but a 
X 2 (likelihood penalty term) was included to tie a given background's branch- 
ing ratio relative to the signal to the measured values within their known 
uncertainties. The yield of D° — > K~fi + i> deduced from a fit to its q 2 -cos8e 
scatterplot served as an estimate of this important background in the fit to 
the q 2 -cos9i scatterplot for D° — > -n~n + v. A more complete description of 
this fitting procedure along with our value for f { +\o)/f{ K) (0) appears m a 
companion paper [9]. 

Figure 4 shows how the q 2 and cos^ projections predicted by the fit compare 
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Fig. 4. The q 2 and cos#£ projections of the two-dimensional fit compared to the 
data histogram for both D° — > ir~ fx + u (upper row) and D° — > K~ p + u (lower row). 
For the case of D° — > 7r~ the peaking backgrounds included the sum of K~fi + v, 
K~tt° /U + i/, K°ir~ p + u and p~ fiu ; while for D° — > K~ fi + v they included K~~tt° p + u. 



to the data as well as the various signal and background components of these 
projections. The results relevant to D° — > K~fi + v are m po i e = 1.93 ± 0.05 ± 
0.03 GeV/c 2 , a = 0.28 ± 0.08 ± 0.07, and f { _ K \o)/f { + K) (0) = -1.7±U ± 0.3. 

The systematic error was determined by comparing results using different 
event selections, alternative fit methods, and looking at the consistency of 
results between split samples. We begin with some of the many alternative 
event selections that were investigated. A f!^\q 2 ) and /^(O) / /^(O) mea- 
surement was obtained from fits where each of these cuts was varied relative to 
our baseline: the detachment of the secondary vertex from the primary vertex 
was varied from 4 to 12 standard deviations, the secondary vertex was required 
to lie out of all target material, the momentum cut on the muon was raised 
from 10 to 25 GeV/c, the secondary isolation cut was tightened from < 1% 
to < 0.1% and the confidence level on the secondary vertex was raised from 
> 1% to > 15%. The split sample compared the form factor information for 
particles to that for antiparticles. Various alternative fits were employed. For 
example, in some fits, the two pole masses were allowed to float while keeping 
f iK \0)/f+ K \0) fixed compared to our standard fit where all three parameters 
were free to float. In another fit variant, the fit was performed on the 5m - q 2 
scatter plot as opposed to the q 2 - cos0£ scatterplot. 

Finally a fit to m po i was made directly from the non-parametric f'^ < \q 2 ) 
results illustrated in Figure 3. This fit minimized a x 2 given by given by Eq. 
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i 3 




where the sum runs over all reporting bins, C 1 is the inverse of the co- 
variance matrix, f + (q 2 )^ are the measured f^iq 2 ) values, and f + (q 2 )^ 
are the predicted f^\q 2 ) within the parameterization. The second term is 
a likelihood penalty term that parameterizes uncertainty in the level of the 
charm background. The parameter b is a background multiplier that multi- 
plies the expected Monte Carlo background yield and cr& is our estimate of 
its uncertainty. The parameterized f + (q 2 )^ depends on a normalization pa- 
rameter /+(0) and a shape parameter m po i c . This fit returned a pole mass 
of m po i e = 1.91 ± 0.04 ± 0.05 GeV/c 2 which is in remarkably good agree- 
ment with m po i e = 1.93 ± 0.05 ± 0.03 GeV/c 2 -the value obtained from the 
two-dimensional, parametric fit. Again, the systematic error of the fit to the 
non-parametric f+\q 2 ) was obtained by checking its stability against a va- 
riety of different fit variants, cut variants, assumed background levels, and 
/-//+ assumptions. 

Using a procedure identical to the parametric procedure used for D° — > 
K-fi+u, we find m po i c value (or D° -> tt>+z/ is m pole = 1.9lt^±0.07 GeV/c 2 . 
The systematic error on this result included an additional important cut vari- 
ant consisting of raising the log-likelihood difference between the kaon and 
pion Cerenkov hypothesis from 3 to 5 and in the process reducing the fraction 
of kaons misidentified as pions by about a factor of two. 



5 Summary 

Figure 5 compares our f+{q 2 ) measurements to a recent [2] Lattice QCD cal- 
culation 6 and our best fit values for m po i c in the pole mass parameterization 
(Eq. 3) and a in the modified pole mass parameterization (Eq. 4). We obtained 
a value of f { _ K) (0)/f{ K) (0) = -1.7±i;| ± 0.3 that is also consistent with the 
value that can be derived from information in [2]. A tabular summary of the 
data of Figure 5 and its correlation matrix has appeared in Table 1 of Section 
3. 

Our fit to the m po i e parameter in pole mass parameterization was m po i c = 
1.93 ± 0.05 ± 0.03 GeV/c 2 . This is compared to previous published data in 

6 We re-scaled their calculations to insure that /+(0) = 1. 
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Fig. 5. The background subtracted /i (g 2 ) (diamonds with error bars) is compared 
to a pole form with m po i e = 1.93 GeV/c 2 (solid curve) , a modified pole form 
with a = 0.28 (dashed curve) , and the unquenched, Lattice QCD, calculations 
given in reference [2] (triangles with no error bars). This form factor is for the 
process D° — » K~ /j, + u. The a and m po i e used for the plots are obtained using the 
two-dimensional, parameterized fit. 



Figure 6. The most recent m po i e is from CLEO [10] who obtain m po i e 
0.05 ± 0.035 GeV/c 2 . All data are remarkably consistent. 
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Fig. 6. Summary of m po i e measurements. All data are consistent with a weighted 
average pole mass of m po i e = 1.91 ± 0.04 GeV/c 2 . The upper solid line shows the 
spectroscopic pole mass at mp*- The lower solid line and two dashed lines represent 
the weighted average and its error. Our weighted average of all data is 5.1 a lower 
than nip* ■ 

Our fit to the a parameter for the modified pole form is a = 0.28 ±0.08 ±0.07 
from the parameterized, two-dimensional fit. This is very consistent with 0.32 
± 0.09 ± 0.07 , the value obtained from our fits to non-parametric data shown 
in Figure 3. The most recent published measurement is from CLEO [10] who 
obtain a = 0.36 ±0.10lo;o7. Our value for the a parameter is 1.9 a lower than 
the value quoted in [2] for D° — > K~fi + u. 7 



We believe that only statistical errors on a are included in [2] 
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We also find that m po i e for D° -> 7r> + z/ is m polc = 1.91±g;fg ± 0.07 GeV/c 2 . 
This value is compatible with our value for the pole mass for D° — > K~fi + u. 
In the naive pole dominance model, the m po i e for D° — > 7T~fi + i> would be at 
the mass of the D* + and would therefore lie lower in mass than rri£,* expected 
for D° -> K-fi+u. 
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